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Abstract 

We consider the Ising model on Z x Z where on each horizontal line {{x,i),x G 
Z}, the interaction is given by a ferromagnetic Kac potential with coupling strength 
J-y{x, y) ^ jJ^j^x — y)) at the mean field critical temperature. We then add a nearest 
neighbor ferromagnetic vertical interaction of strength e and prove that for every e > 0 
the systems exhibits phase transition provided 7 > 0 is small enough. 
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1 Introduction 

We consider an Ising model on the lattice ZxZ, whose points we denote by (x, i). The spins 
a{x, i) take values in {—1, +1} and on each horizontal line, also called layer, {(x, i),x G Z}, 
the interaction is given by a ferromagnetic Kac potential, that is, the interaction between 
the spins at (x, i) and (y, i) is given by 

- y)cr{x, i)a{y, i), ^ v) = 1’ (l-^) 

y^x 

where J^{x,y) = c^'yJ{'y{x — y)); J{r), r G M, is a symmetric probability density with 
continuous derivative and support in [—1,1], 7 > 0 is a scale parameter, is the normal¬ 
ization constant {c^ tends to 1 as 7 —)• 0). We also suppose that J(0) > 0. denotes 
the Hamiltonian with only the interactions dni) on each layer, so that different layers do 
not interact with each other, the system is essentially one dimensional and does not have 
phase transitions. 

We fix the inverse temperature at the mean field critical value /3 = 1 so that also in the 
Lebowitz-Penrose limit no phase transition is present. Purpose of this paper is to study 
what happens if we put a small nearest neighbor vertical interaction 

— e a{x,i)a{x,i + 1). (1-2) 
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The main result in this paper is the following. 


Theorem 1. Given any e > 0, for any 7 > 0 small enough yifj pL~, the plus- 
minus DLR measures defined as the thermodynamic limits of the Gibbs measures with 
plus, respectively minus, boundary conditions. 

It is worth mentioning that a version of Theorem [1] holds for /3 > 1 with e = 7 "^ for any 
A. (See [5] where indeed the above result has been conjectured.) 

In many cases it has been proved that if in the Lebowitz-Penrose limit there is a phase 
transition then in dimension d > 2 there is also a phase transition at small 7 > 0 (i.e. 
without taking the limit 7 —>• 0). We cannot follow this route here because we do not 
know the phase diagram for our model in the limit 7 —)■ 0: a “Lebowitz-Penrose theorem” 
for our system is an interesting open problem that our analysis does not solve. If the 
support of the Kac interaction would contain two dimensional balls (i.e. layers at distance 
of order 7 “^ interact with each other) then the Lebowitz-Penrose analyis [ 8 ] would apply 
and therefore the free energy in the limit 7 —>• 0 would be the convex envelope (i.e. the 
Legendre transform of the Legendre transform) 



where /e(m) is the free energy of the one dimensional Ising model with nearest neighbor 
interaction of strength e. (jl.3p yields a phase transition if e > 0. Does (jl.3p remain valid 
also when the Kac interaction is only horizontal? We do not know the answer but our 
analysis shows that indeed our system has a phase transition as indicated by (jl.3h . 

The proof of Theorem [T] requires a non trivial extension of previous results on Kac po¬ 
tentials and it is given in complete details in this paper. It is obtained by proving Peierls 
bounds for suitably defined contours. The bounds are established via a Lebowitz-Penrose 
coarse graining procedure which however is not straightforward for the reasons explained 
before (due to the local nature of the vertical interaction and the strictly horizontal struc¬ 
ture of the Kac interaction). The trick is to use ferromagnetic inequalities to compare the 
magnetization under with that under the corresponding Gibbs measure for which the 
vertical interaction is removed in a chessboard fashion. To this new system (which is more 
decoupled but not so much as to lose the phase transition) we can apply the Lebowitz- 
Penrose coarse graining strategy. In this way we reduce the proof of the Peierls bound to 
the analysis of variational problems for a suitable free energy functional. 

The model we are considering is related to a d = 1 quantum spin model with transverse 
field, whose hamiltonian is: 

^(<7) = - ^J'){x,y)a^{x)a'^{y) -a^(T^(x) 

x^y X 

in its stochastic representation via Feynman-Kac, 13] and [7j. We are indebted to D. 
Ioffe for pointing out the connection and for useful comments. 

At this point we state two conjectures. 

The hrst question is: can we choose e = 6 ( 7 ) so that 6 ( 7 ) —>■ 0 as 7 —)• 0 and still have for 
all 7 small enough a phase transition? Is there a critical choice for 6 ( 7 )? The conjecture 
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is that setting €( 7 ) = k > 0 , we have a phase transition if k is large enough and no 

phase transition for k small. 

This is related to the next conjecture. Consider the system where on each layer we have 
a process m{r, i) G M, r G M, i E Z. The formal Gibbs measure that we want to study is: 

i 

where P is the Euclidean (f)\ probability measure, namely the stationary solution of the 
real valued stochastic PDE 

du{r,t) = (^-u"{r,t) — u^{r,t)^dt + dw, (1-5) 

dw white noise in space-time. The conjecture is that there is a phase transition for k large 
and no phase transition for k small. 

The measure in (jl.4l) is the formal scaling limit of the Gibbs distribution of the empirical 
magnetization when we scale space as x —>■ r = and renormalize the averages by a 

factor 7 “^/^ as proved in [2] and [ 6 ]; see also [1], where (|1.5I) without the second derivative 
term is derived by studying the critical fluctuations in the mean field version of the model. 

More precisely, in both papers the question is about the analysis of the long space-time 
fluctuations of the d=l Ising model with Glauber dynamics and Kac potential at /3 = 1 
(like ours here). Namely the analysis of the fluctuations held 

X 

with (j) a test function. This is the right normalization because one can prove that at such 
long times the typical values of the empirical magnetization in the limit 7 —>• 0 have order 
7 ^/^ and not the normal values 7 ^/^ of the hnite time huctuations. It is then proved that 
the above huctuations held converges to f (j){r)u{r,t)dr where u solves (II.5p . 

Outline of the paper: As already mentioned, our proof involves the study of the Gibbs 
measures for another Hamiltonian, denoted by and dehned in (j2.3p , where the vertical 
interactions are removed in a convenient chessboard fashion. This yields a two dimensional 
system with pairs of long segments of parallel layers interacting vertically within the pair 
(but not with the outside), plus the horizontal Kac interaction. For this system we can 
exploit the spontaneous magnetization that emerges when two parallel one dimensional 
Kac models at mean held critical temperature interact vertically as in our case, as studied 
in Section [5j This spontaneous magnetization plays a natural role in the dehnition of 
contours (as in the analysis of the one dimensional Kac interactions below the mean 
held critical temperature). The main point is that for this chessboard Hamiltonian, and 
after a proper coarse graining procedure, we are able to implement the Lebowitz-Penrose 
procedure: the corresponding free energy functional is dehned in Section[3]and the problem 
of getting the corresponding Peierls bounds for the weight of contours is transformed in 
variational problems for the free energy functional. This is the content of Theorem [3l 
whose proof involves, as a preliminary step, the study of the free energy function of two 
layers and its minimizers (determining the spontaneous magnetization). There are delicate 
choices of scales so as to allow the implementation of this procedure, as explained in Section 
[2j In Section H] we get an upper bound for the restricted partition function that appears 
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in the numerator of the weight of a contour. Section 0 is dedicated to the study of the 
minimizers of the free energy functional of Section [3l which then will play a crucial role 
in the lower bound for denominator in the weight of a contour, completed in Section [H 
In Section [7] we combine the two estimates to conclude the proof of the main theorem. 
The analysis of the mean field free energy function for two layers and the crucial estimates 
used in Section [5] are carried out in the Appendices. 


2 Contours 


For i S Z, we call i-th layer the set Z x {i}. As mentioned in the introduction we shall 
extensively use coarse graining, for which we start by partitioning each layer into intervals 
of length i G {2", n G Z}. Let , A: G Z} denote the partition of the i-th layer: 


Cf’* = Cix {i} := {[ki, (k + 1)^) n Z) X {i}, where k = \_xli\ (2.1) 


and, as usual, [sj = max{x G Z: x < s}. To simplify notation we restrict 7 to belong to 
{2-^,n G N}. 

For the coarse grained description we shall use three length scales and an accuracy pa¬ 
rameter C > 0 which all depend on 7 : 

7-1/2^ = C = 7“, l>a>a>0. (2.2) 

The smallest scale, 7 ”^'^^, will be used to implement the Lebowitz-Penrose procedure to 
define free energy functionals. Together with (^, the scales I- and £+ will be used to define, 
at the spin level, the plus and minus regions and then the contours. 

For notational simplicity we suppose that also 7 ““ and all the above lengths belong to 
{2”,n G N+}: this is a restriction on a that could be removed by changing “slightly” a 
with 7 . 

We shall prove Theorem [T] for the “chessboard” Hamiltonian; 

= -\'^ J-f{x,y)cr{x,i)a{y,i) - e'^Xi,xCr{x,i)a{x,i + 1), (2.3) 

x^y,i x,i 


where 



if [x/£+J + i is even, 
otherwise. 


(2.4) 


Definition 1. When Xx,i = 1; according to ()2.4h . we say that {x,i) and {x,i + l) interact 
vertically and denote by Vx,i the site (x,j) which interacts vertically with {x,i). 


By the GKS correlation inequalities (see e.g. Theorem 1.21 in Chapter IV of [9]), the plus 
state for is less magnetized than the one for the full Hamiltonian (with Xi,x replaced 
by 1 everywhere). Hence Theorem [T] will follow once we prove that the magnetization in 
the plus state of the Hamiltonian given by (|2.3p is strictly positive. 

For the chessboard Hamiltonian, we shall see via a Lebowitz-Penrose analysis that in the 
limit as 7 —>• 0 there is a spontaneous magnetization equal to some > 0 in the plus 
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state and —m^ in the minus state. This will follow from the analysis in sections [3H51 This 
value is used to define contours, as we now explain (taking mg > 0 for granted). 

Define first the empirical magnetization on a scale .^G{2”,n€N}in the layer i as 


i) 




(2.5) 


We also consider the partition of 7? into rectangles {Q'y{k,j): k,j G Z}, where 
Q^{k,j) = {[k£+, {k + 1)£+) X {j + 1 ) 7 ““)) n Z^ if A: is even 

and 

Q-f{k,j) = {[ki+, {k + 1)1+) X (j 7 "“, (j + 1 ) 7 ““]) n Z^ if A: is odd. 

For convenience we sometimes write Qx,i = Q'y{k,j) if {x,i) G Q'y{k,j). The important 
feature of this definition (frequently exploited in the sequel) is that the spins in each 
rectangle Qx,i do not interact vertically with the spins of the complement, namely recalling 
the definition of Vx,i and that 7 ““ is even, we see that Vx,i G Qx,i for all {x,i). Notice also 
that the rectangles Q'y{k,j) become squares if lengths are measured in interaction length 
units: in fact in such units the horizontal side of a rectangle has length £+j^~^ = 7 ““ and 
the vertical side has also length 7 ““ as the vertical interaction length is equal to 1. The 
other important feature behind the definition of rectangles is that their size in interaction 
length units diverges as 7 —)■ 0 : this will be exploited to prove decay of correlations from 
the boundaries. 

The random variables T]{x,i), 6{x,i) and Q{x,i) are then defined as follows: 

• r]{x, i) = ±1 if \a^^~\x, i) =F m^l < (] r]{x, i) = 0 otherwise. 

• 0(x,i) = 1, [= -1], if ?](y,j) = 1, [= -1], for all (yj) G Qx,i; 0(x,i) = 0 otherwise. 

• e(x,i) = 1, [= -1], if 7(2/, j) = 1, [= -1], for all (yj) G + + 

with (k,j) determined by Qx^i = Q-y{k,j). 

Namely, for the 0 variables we consider a “block” 3 x 3 of Q-rectangles. 

The plus phase is the union of all the rectangles Qx,i such that 0(x,z) = 1, the minus 
phase is where 0(x,i) = —1, in the complement the phase is undetermined. 

Two rectangles Q^{k,j) and Q'y{k',j') are said to be connected if {k,j) and {k',j') are 
♦-connected, i.e. |A: — A:'| V \j — j'\ < 1. By choosing suitable boundary conditions, we shall 
restrict in the sequel to spin configurations such that 0 = 1 outside of a compact (the 
case when 0 = — 1 can be recovered via spin flip). Given such a a, we call contours the 
pairs T = (sp(r), ? 7 r)) where sp(r) is a maximal connected component of the undetermined 
region, called the spatial support 0 /T, and rjr is the restriction of 7 to sp(r), called the 
specification o/T. 

Geometry of contours. Denote by ext(r) the maximal unbounded connected compo¬ 
nent of the complement of sp(r) and 9out(r) the union of the rectangles in ext(r) which 


5 


are connected to sp(r). 9in(r) is instead the union of the rectangles in sp(r) which are 
connected to ext(r). 0 is constant and different from 0 on (9out(r) and we call plus a 
contour T when 0 = 1 on 9out(r) and minus otherwise. Observe that in a plus contour 
T] = 1 on 9in(r). 

Analogously we call intfc(r),A: = the bounded maximal connected components 

(if any) of the complement of sp(r), (9in,fc(r) the union of all rectangles in sp(r) which 
are connected to intfc(r). 9out,fc(r) is the union of all the rectangles in intfc(r) which are 
connected to sp(r). Then 0 is constant and different from 0 on each 9out,fe(r) and we 
write (9^^ (^) fe(^) if © = ±1 on the former, observing that 77 = ±1 on 

^ (r) , respectively. We also call 

c(r) = sp(r) u [Jintfc(r). (2.6) 

k 


Diluted partition functions. Let A be a bounded region which is an union of Q- 
rectangles. The plus diluted partition function in A with boundary conditions a is 

=: ^A,^(0 = Ion din(A)), (2.7) 

O-A 

where d is a configuration on the complement of A; 0 is computed on the configuration 
(o'Ajd) and din(A) is the union of all Q-rectangles in A connected to A'^. Minus diluted 
partition functions are defined analogously. As a rule we denote by Z\^a{A) the partition 
function with the constraint A, A & set of configurations. Notice that there is no vertical 
interaction between the spins in A and those in its complement because A is union of 
rectangles. 

The plus diluted Gibbs measure (with boundary conditions a) is defined in the usual way, 
namely, given a configuration of spins ua on A, the weight assigned to ua by the plus 
Gibbs measure is given by 

^A,, 7 (^a) = --l{0=lon Sin (A)}- (2-8) 

^A,ct 

The minus diluted Gibbs measure is defined analogously. 

We shall prove the Peierls estimates for the plus and minus diluted Gibbs measures, which, 
as a consequence, have distinct thermodynamic limits; Theorem [1] will then follow. 


Weight of a contour. We are now ready to define the fundamental notion of weight 
of a contour. Let T be a plus contour (the definition for minus contours is obtained by 
spin flip) and a a configuration on the complement of c(r) such that r/ = 1 on (9out(r) (in 
agreement with the definition of a plus contour). Then the weight of T with boundary 
conditions a is 


Wr{a) 


Zc{ry,a{v = Vr on sp(r); 0 = ±1 on each 

Zc{ry,a{® = 1 on sp(r) and on each5^^^ ]^(r)}) 


(2.9) 


where Z\^^{A) is the partition function in A with Hamiltonian with boundary con¬ 

ditions a and constraint A. In the next sections we shall prove the following theorem 
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Theorem 2 (The Peierls bounds). There are c > 0, eg > 0 and 7 . : (0,oo) —>■ (0, 00 ) so 
that for any 0 <e<eo, 0 < j < and any contour F with boundary spins d 

Wr(d) < (2.10) 

In Section [7] we shall see how to prove Theorem [1] using the Peierls bounds (j2.10l) . 


3 Reduction to a variational problem 


The goal of this section is to introduce the Lebowitz-Penrose free energy functional and 
to set the variational problem that emerges in the estimates of the partition functions in 
(EJl). We start by the next proposition which deals with the very simple situation of two 
layers of ±1 spins whose unique interaction is the nearest neighbor vertical one. It is just 
a chain of independent pairs of spins. Therefore the multi-canonical partition function, 
where we fix the magnetization on each layer, is studied by very simple tools. This first 
result, proved in Appendix lAl for sake of completeness, describes its convergence (in the 
thermodynamic limit) to the infinite volume free energy 4 >^(mi,m 2 ) and finite volume 
corrections. We then state and prove a Lebowitz-Penrose theorem for the spin model 
associated to the chessboard Hamiltonian 


Proposition 1. Let n be a positive integer, Xn = { — 1,1}""- For i = 1,2, let 
mi G {—1 ^ : j = 1 ,... , n — 1 } and set 


Ze,n(mi,m2) 


E 






(3.1) 


There is a continuous and convex function defined on [—1,1] x [—1,1], with bounded 
derivatives on each [—r,r] x [—r,r] for |r| < 1, and a constant c > 0 so that 


~ log Tl 1 

- 4 (mi,m 2 ) - c- < - logZe,„(mi,m 2 ) < -(/>e(mi, m 2 ). 

n n 


(3.2) 


We shall next use the above proposition to study the partition functions which enter in 
the definition of contours. We thus consider a region A which in the applications will be 
the spatial support of a contour. Here it only matters that A is a connected set union of 
Q-rectangles. We want to bound from above and below the partition function 

Za 7 (^) := ^ (3.3) 

(ja&A 

where d is a spin configuration in the complement of A and “the constraint” .A is a set of 
configurations in A defined in terms of the values of r]\. We shall coarse-grain on the scale 
\Ye thus call M^- 1/2 the possible values of the magnetization densities 
has been defined in (12.511 ). namely 

= {-1, -I + 27 I/ 2 ,..., 1 - 27^2, 1 } 

and we set 

WIa := {m(-) G (M^- 1 / 2 )^ : m(-) is constant on each C'^ ^ C A}. (3.4) 
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The Lebowitz-Penrose free energy functional (on A with boundary conditions m) is the 
following functional on [— 1 , 1 ]^ (whose elements are denoted in short by m) 

1 




(a:,i)GA 


- ^ J^{x,y)m{x,i)m{y,i) 

J'y{x,y)m{x,i)fh{y,i), 

{x,i)eA, {y,i)ik 


(3.5) 


where fh G [—1,1]^'^, </>e is the free energy function in (|3.2p and Vx^i is given in Definition 
[TJ (Recall that Vx,i G A for each (x, i) G A since there are no vertical interactions between 
a Q-rectangle and the outside.) 

Notational remark. The same formula is used when m is defined in a set A contained in 
the complement of A; in such a case the sum over {y, i) in the last term is extended only 
to A. 


By an abuse of notation we write, analogously to (12.5p . 


(3.6) 


and define r]{x, i; m) = ±1 if |m*'^“^(x, z) =F < C and = 0 otherwise. We still denote by 
A a constraint that depends on r]{-] m) as for instance m) = on A. 


Theorem 3. There is a constant c so that 


\ogZx{<7]A) < - inf TA_.y(m|m) + c|A| 7 ^/^ logy \ 




(3.7) 


_ 1/2 

where, recalling (12.51) . fh{x,i) = a'* {x,i), {x,i) ^ A. Moreover, for any m G A4a n A 


logZA(<T;^) > -FA, 7 (m|m) - c|A| 7 ^/^ logy ^ 


(3.8) 


Proof. Here is essential the restriction to regions with no vertical interaction with the 
complement. We have, writing a for a spin configuration in A, 

ZA{a;A)=Y, E 

(■)=m(.) on A 

where is the Hamiltonian ()2.3I) in the region A interacting with a outside A. 

By the smoothness of we get 

I 5 E E(a(x ,y) - J^{x,y)^cr{x,i)a{y,i)\ < c|A| 7 ^/^ 

{x,i)eAy^x 

where cr{y,i) = a{y,i) if {y,i) ^ A and 

1 


j.^{x,y)-.= ^ E E 
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x'&c2 


Thus, recalling that there is no vertical interaction between A and its complement, 




{x,i)eA 


where is the Hamiltonian with e = 0 and the spins replaced by 

m{x,i) and fh{x,i). We then get, using (|3.2|) . 


ZA(a;A) < IA^aI sup e 2 l^(^,i)eA 

meA 


<pe(m(x,i),m(v:^^i)) cjAjA^^ 
^ 1 


which proves (I3.7p because |A4 a| < (cy (for a suitable constant c). (13.81) is 

proved similarly. □ 


The variations of on the scale 7 give a contribution of the order | A| 7 ^/^ to the errors 
in dsm and (|3.8p : in (|3.7I) there is also a contribution of order |A| 7 ^/^ log 7 ^ coming from 
the cardinality of A4 a- In (13.8p we need to take into account the lower bonnd in Proposition 
[TJ Of conrse in the upper bound of the partition fnnction we can drop the condition that 
m takes values in M^- 1/2 and that it is constant in the intervals C'^ ’* C A: 

Corollary 1. In the same context of Theorem\^ 

\ogZ\{a;A)<— inf Fa 7 (ni|m) + c| Ajy^/^ log 7 “^. (3.9) 


4 The upper bound 

Let us now be more specific and see how ()3.9p is used to get an upper bound for the 
nnmerator of ()2.9I) . The key point will be to prove that the excess free energy due to the 
constraint 7 = r/r is much larger than the errors in (I3.7I) - (I3.8I) . 

In the sequel we specify A = sp(r), and refer to the paragraph “Geometry of contours” 
in Section [2j Because the notation gets clnmsy in some formnlae, we shorten it a bit as 
follows: 

Ain = 5in(r), = int^(r), 

recalling that the suffix ± here refers to the (constant) valne of 0 on the corresponding 
9out,fc(r), and that 7 = ±1 on A^. Set then 

Ao = sp(r) \ (Ain U {UfcA+} U {UfcA^}) 

so that one has the following partition of c(r): 

c(r) = Ao U Ain u {UfcA+} U {UfcA-} U {Ufc/+} U {Ukl^}. (4.1) 

Thns the fnnction fh in (13.9p is specified by the spins outside c(r) and by those in the sets 
When necessary we write tha ^ for its restriction to the complement of c(r) 

^ k 

and to respectively. Finally the constraint in Theorem [3] is A = {i] = rjr on A}. 
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Following m (see Chapters 6 and 9), an important ingredient in the proof of Peierls 
bounds consists in showing that the minimizers of the free energy functional have good 
regularity properties even when constrained to have given magnetization values in small 
boxes (the multi-canonical constraints). The next proposition, proved in Appendix ICl 
shows that the infimum in (l.f-hh can be restricted to smooth functions. 


Proposition 2. There is a positive constant c so that, with the same notation as above 
and recalling that A = sp(r), 


inf F/^y{m\m) = inf (4-2) 

where 

S/\Q\={m\ sup \m{x,i) — m^^~\x,i)\ < c^°^}. (4-3) 

(a;,i)eAo 


Remark. The smoothness request could be extended to the whole A without changing the 
infimum but we only need it in Aq. 

Let us write 


— m{x, i)m{y, ~ 2 ( ~ + [m{x, i) — m{y, 

in some of the terms. 

With the above notation, and recalling ()3.5I) . we get 

Fsp{r),'yim\fh) = T^o,.^(mAo) + T^i„,^(rnAiJm^ext) 

+ E ^ 1 +) + E 1 ^-,- ) 


k 


+ h E E -m{y,i)f, 


(x,i)€Ao {y,i)^Ao 


where, writing m for 


^Ao.7M = E + 

(x,i)GAo 

+ 2 J-i{x,y){m{x,i)-m{y,i)f, 




while writing m for ttia- 


-P’Ai„,7("^l"i<^ext) = ^Ai„7("l|m^ext) - E 


m[x, i 


{x,i)eAi, 


where 

(^x,i ^ ^ J'yiFiV')- 

y:{y,i)€Ao 

F' ± (m^±\ma ±) is defined analogously. 

^k k If. 


(4.4) 


(4.5) 


(4.6) 


(4.7) 


(4.8) 
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Remark. Recalling the remark after (13.51) . notice that the interaction between Ajn and 
Aq present in Fgp(r), 7 (?TT-|th) is absent in It is instead contained (and 

split via (j4.4l) i in the following terms: (i) the last term on the right hand side of (I4.5p . 
(ii) the second term in the right hand side of (|4.7p with {x,i) G Aq. Something 

analogous takes place to the interaction between A^ and Aq, involving the last term on 
the rhs of (4.4) and the corresponding terms in F'^± .^(^a=*= )■ 

By Proposition [2] and (|4.5I) . we get, dropping the last term in (14.51) . 

Corollary 2. 


where 


inf Fsp(r)y{m\m) 


> ^Ao + ^ Ai„ ) + Y1 ^ A+ i^^r+ ) 

k * 

k * 


(4.9) 


0 

0 

= inf < 

F^o.7("i) ?4iG[-l,l] 

^in 

(™lText 

= inf 

^Ai„7("l|™lText 

+ 0 

ilFa^+] 

^k 

= inf < 


t> 1 

1 

{rria ) = inf ■ 



, \m — ^1 < C 7 “, r]{-; m) = 

m G [-l,l]^'",77(-,m) = 1,}, 
m G , 7 (-,m) = l|, 

m G [—1, , T/(-, m) = —1|- 


}. 


(4.10) 


Corollary [2] is useful for us because it allows to split the original variational problem on 
the left hand side of (|4.9I) into separated, localized variational problems, as on the right 
hand side of (j4.9h . 

Recalling (13.91) . we have the following upper bound for the partition function in the nu¬ 
merator of ( 1 ^ : 

g-$Ao+dA|7l/Rog7-lg--I-Ai„(m.,,t){JJ Z+(/+)}{[] Z-(/-)}, (4.11) 


where 


z+(4P = E' 


-H(» +)-»7 
h A 




(4.12) 


The superscript -|- in the sum means that the sum is restricted to spin configurations in 
it such that a configuration made by aj+ in it and by any configuration with 7 = 1 in 

A+ has 0 = 1 on 5+^^;. (P), see (12.9p . Z {ly. ) is defined analogously. 

Following the Peierls strategy we use at this point the spin flip symmetry to rewrite ()4.1ip 
in a more convenient way. In fact we have: 


) = 4>+ (m_^ ) 

^k -'fe ^k ^k 


(4.13) 


and therefore Z (7^. ) = ). The numerator in (12.9p is thus bounded by 

= IT on sp(r); 0 = ±1 on each 

< g-$Ao+c|A|7i/2log7-i 




{n244P){n^44-)}. 


(4.14) 
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The key point is now to prove a lower bound on <l>^p so good as to kill the error terms in 
the first exponent and to give what is required by the Peierls bounds. The other factors 
in (I4.14P will simplify with those coming from the lower bound modulo a small error. 
Preliminary to that is the analysis of the two layers free energy f^{mi,m 2 ). In Appendix 
m it is proved that: 


Proposition 3. For any e > 0 small enough 


fe{mi,m2) := 

has two minimizers, := and there is a constant c so that 

\me — V^\ < 

Moreover, calling /^^eq the minimum of fe{m), for any ( > 0 small enough: 

fe{m) - A,eq > , for all m^U^, 


where 


(4.15) 


(4.16) 


(4.17) 


Ur := 


m 2 ) : \mi - me\ < ^, i = 1,2| U |(mi,m2) ; |mi + me| < |, i = 1,2|. 


C . 


The following lower bound for <I>Ao follows from Proposition [3] and proves that the excess 
free energy <I>Ao — /e,eq^^ grows at least like c|(recall that f^^eq is defined after 
(j4.16p and that = /e,eq-^^ when m is identically equal to rUe or to —mg). As 

desired such excess free energy c|Ao| 7 '^“^^“ is much larger (for small 7 ) than the error 
term in the first exponent in (|4.14l) which is given by c|sp(r)| 7 ^/^ log 7 “^. 


Theorem 4. There is c> 0 so that 

3>A,> /..=,^ + C 1^" 




7-(^-“) min{7";72“}. 


(4.18) 


Proof. We rewrite (j4.6p as 




(x,i)eAo 


+ i ^ J^{x,y){m{x,i) - m{y,i)^ (4.19) 


and start by bounding from below the first term. We distinguish two cases: 


• (i) when ? 7 (x,z) = ii{vx,i) / 0 we bound fe{m{x,i),m{vx,i)) > A,eq; 

• (ii) in the other cases for all {y,i) G {m{y,i),m{vx^i)) ^ for 7 small 

enough (as, by the smoothness condition, |m — < 07 "). We then bound 

/g(m(x, i),m{vx, i)) P fe,eq + cC^. 
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Thus the hrst term on the right hand side of (|4.19l) is bounded from below by 


/e,eq^^ + iV(ii)C7 

where, writing Vx^i = is the number of distinct pairs of intervals ci~’^ 

where case (ii) occurs. 

The second term on the right hand side of (j4.19p is bounded from below by retaining only 
the terms where {x,i) and {y,i) are in two consecutive C^~’^ intervals and r]{x,i;m) = 
—rj{y,i;m) / 0. Suppose for instance rj{x,i]m) = 1, then \m{x,i;m) — me| < 2Q (for 7 
small enough and using smoothness as before). Analogously \m{y, i) + m^\ < 2C, and recall¬ 
ing the assumption that J(-) is strictly positive at the origin, we see that the contribution 
to (I4.19P coming from any such pair of intervals is, for 7 small enough, at least 

C7^(7“^^““^)^ = c7-(^-")7“. 

To conclude we observe that by definition of contours for any Q-rectangle Q' in Aq there is 
a rectangle Q in Aq connected to Q' (or just Q' itself) with the following property. Either 
case (ii) occurs in Q or there are two consecutive intervals one at least inside Q with 
opposite values of y, or both the above events occur. (|4.18l) is then obtained because the 

lAol 

number of Q rectangles in An is — 7——r-. □ 


5 Characterization of minimizers 

The lower bound of the denominator in (12.9p will be obtained by computing the free energy 
functional on a suitable test function m on sp(r). 

On Aq, m will be a constant approximately equal to mg while it will be (approximately) 
equal to the minimizers of ^Ain(tho-ext)) +)) ^\-{rh-a _), in the respective 

sets Ain, A^, and A^. 

The main difficulty will be to estimate the last term in ()4.5p which in the upper bound 
for the partition function could be neglected being non negative. We shall prove below 
that the term [m{x,i) —m{y,i)]^ in (14.51) with (x,i) G Aq and (y, i) ^ Aq is bounded from 
above by e~^'^ , with the above choices of m{x, i) and m(y, i). Thus the last term in (14.51) 

will then be negligible also in the lower bound. 

The analysis of the minimizers is essentially the same for all of them and, for the sake 
of definiteness, we will just look at the minimizer of <1>J. (m-o-ext)) referring to (jd.lOp and 
(j4.7l) for the definition. 

Recalling (14.7p . (13.5p and (14.8p we have 

^(j)e{rn{x,i),m{vx,i)) - ax,irn{x,i)'^'j 

{x,i)eAin 

- ^ J-tix,y){m{y,i) + ^ J.^(x,y)m(y,i)^ (5.1) 

(a;,i)eAin (j/,i)GAin {y,i)tsp{r) 
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Being a continuous function of {m(x, z), (x, i) G Ajn} the function m i—)• F'^, ^{rn\m) 
attains a minimum when m varies in the compact set 

K := Pi — me| < c|- (5.2) 

(a;,i)GAin 

We are going to prove that the minimizer is unique and will establish properties of the 
minimizer typical of the correlations in the Gibbsian high temperatures regime. 

We fix arbitrarily a pair (a:,i) and = Vx,i of vertically interacting sites in Ain, and 

regard ^{rn\m) in (15.ip as a function of m{x,i) and m{x,i') alone with all the other 
m{y,j) considered as fixed parameters that we denote by u{y,j). Let 

= \^u{y, j),y ^ x,j = i,i': {y,j) e A'^, u{y,j) e {-1,1) and 

\u^~{yJ) - m-el < C, for all {y,j) i (5.3) 

namely the set where the function ?/(-;zi) is identically 1 except maybe on the intervals 
containing {x,i) or {x,i') where we do not impose conditions on the tt(-). 

For any u G we introduce the function 

ge{mi,mi>) := 4(m-i,mi') - ^{aimf + ai>ml) - - A“mj/ (5.4) 

where {mi,mii) G (—1,1) x (—1,1), aj is a shorthand for axj, j = and 

Xj = ^ J^{x,y)u{y,j), j = i,i'. (5.5) 

In Appendix [D] we shall prove: 


Proposition 4. There are cq > 0 and 7 . : (0, 00 ) —)• (0, 00 ) such that for any 0 < e < eo 
there are r < 1 and coefficients Cx,i,ii{j,j'), j,j' G so that the following holds for 

all 7 < 7 e. 


For any u G there is a unique minimizer of g^- 

^ Cx,i,i'{j,j') < r for j = i, i'. 


I A“ — A^? I 

< Cx i i'{j,j ')— - - - —, j = i, i', where A®^ is the value of \j when 

’ ’ 1 — a ■' •' 


I (■“) 

\m- — me| ^ 

u is identically equal to m^. 


•x,j' 


|A« - Afl 


1 — a 


< C + C 7 " for j = i, i' 




For any u,v £ Afx i v and j = i, i', 


my _ 


\\U _ \v I 

— / ^ i J ) 
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As a consequence of the above proposition we have: 


Theorem 5. In the same context of Proposition and for 7 small enough the following 
holds. There is a unique minimizer m* of F'^, in K, see (1521) and for any {x, i) G 

Ain 

\m*{x,i) — m^l < C, (5.6) 

\m*{x, i) — < 2r'^, (5.7) 

where n is the minimal number of steps required to go from (x, i) to the complement of 
sp(r) when horizontal steps have length < 7 “^ while the vertical steps have length 1. 


Proof. We shall preliminary prove that for 7 small enough the minimizer in Proposition 
IHsatishes — m^l < (. Indeed: 


m 


(u) 


— me 




|a: 


^ 

^ ■ < rC + C 7 " < C, 


1 — a 


(5.8) 


having used the bounds on and |A“, — A^l stated in Proposition [H The last 

inequality + 07 “ < C holds for 7 small enough, because r < 1 and by the choice of f 
and a. 

Since ^{m\fh) is a continuous function of the coordinates m{x,i), {x,i) G Ain, h has 

a minimum in the compact set K. Let m be a minimizer, and ci~’^ a segment in Ain 
whose points are denoted (xi, i),.., (xat, i). Let mx^-i be the function obtained from m 
after replacing the elements m(xi,i) and m{vxi^i) by the minimizer of relative to the 
points (xi,i) and Vxi,i and with u = m on the complement of {{xi,i),Vx^,i\. We then 
dehne iteratively the sequence mx^^...^xk-,i^ k < N, hy applying the above procedure to 
AVe claim that mxi,...,XN\i — 'yi.''^xi,...,xk;i\lh) is non increasing 

in k (because we are relaxing the condition 7 = 1 in and because we are putting at 
each step the minimizer of the corresponding 7 ^) and therefore 


F'Ai^,'yi'^xu...,XN;i\m) < F^.^_^(m|m). 

By ([52]) \mxi,...,XN-AyJ)-'^<i\ < C for all {y,j) in df"/*U. Therefore mxi,...,XN;i ^ K 
and since m is a minimizer 




(m|m) < F'^ Jm., 




m] 


This means that at each step 

FA^,A'^xi,...,xk-cAkh) = FA^^^{mx^,...,xk-Akh) 

and by the uniqueness of the minimizer of g^ we conclude the proof of the claim. Observe 
that we have also proved a sort of DLR property, namely that if m is a minimizer then its 
values at (x,i) and (x,z') minimize the corresponding 7 ^. 
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By the arbitrariness of the choice of (x, i) G Ai^ in the above argument we deduce that 
for all (x,i) € Ain \m{x,i) — me| < C, and since {m{x, i), m{vx,i)) is the minimizer of the 
corresponding then by the second property in Proposition [H and by (15. 5 p 


m{x, i) — rrie 




1 




- md, 


(5.9) 


where m = m outside sp(r). The inequality (15.91) can be iterated n times with n as in the 
text of the theorem and we get 

\m{x,i)-m^\ < X] Cx,i,i'{i,ji) X^ 

yiy^x:{y,j)^Ao 

J. QjlJ-i jiy 

j2=ji,j[ yii=yv-{yi,h)i^a 


We bound the last factor by 2, the sum over the yk is normalized to 1 hence |m(x, i)—me\ < 
2 r". 

It remains to prove the uniqueness of the minimizer of F'^, in K. Suppose there 

are two minimizers m and m', then by the last statement of Proposition [4] 


\m{x,i)-m'{x,i)\< Cx,i,i'{i,j') YZ \m{y,j') - m'{y,j')\. (5.10) 




^X,J' 


The inequality can be iterated n times. But now n is arbitrary because m{y,j) = 'rn'{y,j) 
outside sp(r) and therefore m{x,i) = m'{x,i). 

□ 


6 The lower bound 


In this section we will prove a lower bound for the denominator in ()2.9I) . 

We call trial function a function m € A4sp(r), see (13.4p . namely with values in M^- 1/2 and 

constant on the intervals ’* contained in sp(r). Denote by a the collection of spins 

in c(r)'^ and in the sets 1^ For any such a we choose a trial function mo- and using (13.81) 
we get that the denominator in (12.9p is bounded from below by 


g-c|sp(r)|7i/2- 


/ ^ e k k 


( 6 . 1 ) 


having used the notation of (|4.12p . The whole point is now to reduce (|6.ip to (|4.14|) and 
this is done with a good choice of the trial function. We fix a and start with the function 
m' which in Aq is identically equal to m^, in Ain h is the minimizer of <I>J, (m-o-ext) ™ 
it coincides with the minimizer of <I>^7(?7io-). Since the values of m' are not necessarily 

k 

in M^- 1/2 m! may not be a trial function. We then define m" which at each (x, i) is equal 
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to a value in M^- 1/2 which minimizes the distance of m'{x,i) from M^- 1 / 2 . m" may not 
be constant on the intervals C'^ so that we define rrifj as 

^ m"{y,i). (6.2) 

rua is a trial function and we can use it in (16.Ij) . We claim that 

Fsp(T)ai'<^<T\fha) < Fsp(r),^(m"|?fi<^) + c7^/^|sp(r)|. (6.3) 

Proof. Recalling the definition (13.Sp of the free energy functional, we observe that by 
convexity 

^e{m”{x,i),Tn"{vx,i)) > ^ ^^{ma{x,i),m^{vx,i)). 

(x,i)Gsp(r) (a;,i)esp(r) 


For the terms in (j3.5l) that contain Jy, replacing rria by m" gives an error that is bounded 
from above by c 7 ^/^|sp(r)|, yielding (j6.3p . Similarly 


Fsp{r),-f{m''\ma) < Fsp{v)n{^'\^<^) + C 7 ^'^^|sp(r)|. 

So far we have proved that the denominator in ()2.9I) is bounded from below by 

-Rp(r) ,7 ("*'|m.a) - Efc(0-^+) - Efc(o-^-) 


,-c|sp(r )|772 


E 


(6.4) 


(6.5) 


k k 


(with c a suitable constant which takes care of all the above errors). We next use (14.5p to 
get 

Fs^{r)ai^'\ma) < 

where we used that (i) m' = in Aq; (ii) m' is the minimizer of 4>^±(mo-) and of 
<1>^. (rha) in the respective sets; (iii) the last term in (14.5|) is bounded using (j5.7p . 

In conclusion 


^c{r)-,a{v = 1 on sp(r); 0 = ±1 on each cl^(r)) 

> g-/e,eq^-c(|sp(r)|7l/2 

X Z+(/+)}{J] Z+(4-)}. (6.7) 


Proof of Theorem [B A comparison with the upper bound (I4.14p and use of (I4.18P then 
completes the proof of the Peierls bounds. □ 
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7 Proof of Theorem [T] 


The proof of Theorem [T] is based on the validity of the Peierls bounds (12.101) and it follows 
closely the well known proof for the nearest neighbor Ising model at low temperatures. 

Let {A„} be an increasing sequence of bounded Q-measurable regions which invades the 
whole space and - plus diluted Gibbs measures with boundary conditions a (a may 
depend on re). We want to prove that for 7 small enough and all boundary conditions a 
as in the paragraph of ( 12 .9p 


lim ret - 


0 ( 0 ) < 1 


1 


(7.1) 


By the definition of plus diluted Gibbs measures the event in (17.ip can only occur if there 
is a contour T such that the origin belongs to c(r). Gall A^(r) the number of Q-rectangles 
contained in sp(r). Then there is a horizontal translate of sp(r) by k£^, k < A^(r), so 
that the translate of sp(r) contains the origin. 

This means that 




0 ( 0 ) < 1 




u u {Tfc is a contour} 

r:sp(r)30fc<Ar(r) 


(7.2) 


where T^ is the contour obtained from T by translating it by ki~^. By subadditivity, the 
right hand side of ()7.2p is bounded above by 


r:sp(r)90fc<Af(r) 


Tfc is a contour 


(7.3) 


Now the probability inside the double sum in (|7.3p equals 
^c(rfc);s(?? = Wk on sp(rfc );0 = ±1 on each 






(7.4) 


the inequality justihed by the fact that the denominator in (|2.9p is bounded above Z~^ - 
(since the sum defining the latter quantity contains the terms in the sum defining the 
former one). 

Therefore by (I2.10p . and using the fact that |sp(rfc)| = |sp(r)|, we find that 




0 ( 0 ) < 1 


r:sp(r)90 


(7.5) 


On the other hand |sp(r)| = A^(r )7 " so that the sum on the right hand side of 

(j7.5p is just the sum over all connected regions D made of unit cubes of 



0 ( 0 ) < 1 


< 


E 

D90 


L»|e 


—c|r)|7“i+2“+2“ 


(7.6) 


Since a and a are much smaller than 1 , then the sum vanishes in the limit when 7 —)■ 0, 
see for instance Lemma 3.1.2.4 in uni, so that dni) is proved. 

By (17.11) and the spin flip symmetry it follows that there are at least two DLR measures, 
hence by ferromagnetic inequalities the plus and minus DLR measures of Theorem [T] 
are distinct and Theorem [1] is proved. There are many more consequences of the Peierls 
bounds, see for instance Ghapter 12 in m, but we shall not discuss such extensions here. 
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A Proof of Proposition [T] 


We prove Proposition [T] via equivalence of ensembles. The grand-canonical partition func¬ 
tion TTg of the two layers is trivially equal to the logarithm of 

_ 2 |ge cosh(/i+) -|- e~^ cosh(h_)}, 

(o-i,cr2)G{-l,-|-l}2 


where = hi + and /i_ = /12 — hi. Thus the pressure is given by 

7 re(/i+,/i_) = log(2Z), Z = {e*^ cosh(/i+) -b cosh(/i_)}. 


(A.l) 


We can easily check that the function (h+,/i_) 1 —)• 7r£(/i+,h_) is strictly convex, namely 
its Hessian, denoted here by is a positive definite operator. Indeed, by computation 

we have: 

d . ,sinh(/i_|_) d _gSinh(h_) 


dh. 


-7re(/i_|_, /i_) = 


dh- 


Z 


(A.2) 


9^ 

= z 

52 


(e«h(M-«- 2 ^)> 0 . 


9h2 

52 


I- — e 


6 / 

7r,(/i+,h_) = — (^cosh(/i_) 

sinh(/i+) sinh(h 


gSinh 2 (/i) 


Z 


) > 0 , 


TTe{h+,h-) = -- 


dh+dh-“"^'"^'‘"~' Z 2 

It then follows that the diagonal elements of D‘^'K^{hj^,h-) and its determinant, given by 
|T> 27 re(/i+, h-)\ = + 2cosh(2e) cosh(h+) cosh(/i_) -|- cosh2(h+) cosh2(h_) 


— sinh2(h_|_) sinh2(h_)^. 


(A.3) 


are all positive, and therefore the 2x2 Hessian matrix is positive definite. We now consider 
the Legendre transform of vr^: 


where m = and 


and we have: 


(j)eim) := sup (]-{h,m) - vr, 

h ^2 


(h, m) := h+m+ + h-m-, 


(A.4) 

(A.5) 


Lemma 1. For any m = such that \mi\ < 1, i = 1,2, where 

mi = (m_|_ — m_)/2, m 2 = (m_|_ -|- m_)/2 


(A.6) 


the sup in (IA.4h is a maximum, achieved at a unique h = (/i+,/i_), which is the unique 
solution of 

d d 

m+= 2-—7re(/i+,h_), m_ = 2-—7r,(h+,/i_). (A.7) 


9h_i 


dh- 


In other words, for this h 

4>e{m) = - 7re(h) = sup {^{h' ,m) - 'Ke{h')^. 


(A.8) 


19 
















Proof. If |mj| < 1, i = 1, 2, the function 

T{h) := ^{h,m) - TT^{h) 

goes to —oo when \h\ — )• oo. Together with the continuity of r(/i) this implies that the 
supremum is a maximum, achieved at the critical point, hence ()A.7p . Uniqueness follows 
from the strict convexity of tt^. 

□ 

The following lemma is an immediate consequence of the strict convexity of and the 
properties of the Legendre transform. 

Lemma 2. is strictly convex. Writing D(j)^ for its gradient, m = solves the 

Q 

equation D(l)^{m) ~ if’ ^ ~ */ only if 

m = 2Dtt,.{9). (A.9) 


More explicitly: 


m+ = 2e^ 


sinh(0+) 


cosh(0+) + e cosh(0_) 


171- = 2e 


sinh(0) 


e*^ cosh(0+) + e~'^ cosh(0_)' 


(A.IO) 


Changing back to coordinates (/ii,/i 2 ) and (mi, m 2 ), let ^e(/ii,/i 2 ) and 0e(mi,m2) denote 
the functions 7re(/i+,/i_) and (/>e(m+,m_) when h± and m± are expressed in terms of 
(/ii,/i 2 ) and respectively (mi,m 2 ). Thus vTe and (/>e are the Legendre transform of each 
other: 

(^e(mi,m 2 )= sup ((/I'lmi+ /i 2 m 2 )-:^e(/ii,/i 2 )) = (/iimi+ /i 2 m 2 )-vre(/ii,/i 2 ) (A.11) 

{KK) ^ ^ 


where in the last equality {hi,h 2 ) are functions of {mi,m 2 ) via ()A.7p . 

The following lemma is an immediate consequence of the above. 

Lemma 3. The map {mi,m 2 ) (/>e(mi,m 2 ) is strictly convex on (—1,1)^ and m = 

{mi, m 2 ) solves the equation D(j)e{'hh) = 0,6 = ( 01 , 02 ), if and only if 

m = DTT^{6), (A.12) 

where Dire and Dcpe denote the gradient of ire and 4>e. Moreover, calling 

m = (m+,m_), m+ =---, m_ =---, 


0 = ( 0 +, 0 _), 0 + 


01 + 02 


we have that D(f)e{m) 


if and only if D(fe{m) 


0 

2 ' 


O2-O1 

2 
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Proof of Proposition [T] 

It remains to prove that (/>e is the canonical free energy of the two layers system. 
With (hi, / 12 ) as in (lA.lip and Z as in (lA.ip . we get 


= e = TTijn, z = 1,2 , 

X=\ 


(A.13) 


where Phi,h 2 ,e is the Gibbs measure with Hamiltonian 

n n 

-^{hi(Ti(3:) + h2(T2{x)} - e'^ai{x)a2{x) 

X=1 X=1 

hence the upper bound in (|3.2p . The lower bound follows by an application of Lemma|3]be- 
low, which gives a(n elementary) local limit theorem lower bound for the product measure 
Phi,h 2 ,e of the form const with a uniform constant over n and m. □ 


Lemma 4. With notation as above (see IIA.13\) ). there is a constant c > 0 independent of 
m = {mi,m 2 ) G { — 1 + , 1 — and n> 1 such that 


P, 


h\,h2,e 


ai{x) = miU, i = 1, 2 


.X = l 


> cn 


-zji 


(A.14) 


Remark The sharp bound is cn 
shorter argument, given next. 


but iA. 14 \) is enough for our purposes, and requires a 


Proof of Lemma We first note that '^x=ii^ii^)W 2 {x)) is a two dimensional ran¬ 
dom walk with jumps to the nearest diagonals. We rotate by —vr/d radians and rescale 
space by l/\/2 in order to get a simple 2-d random walk, denoted by A" = {X{n) = 
(Xi(n), X 2 (n)))n> 0 ) with mean jump m* = (m+,m_)/2. Notice two things: 


1. m* G {{i,j)/n : i,j = —n + 2, —n -|- 3,..., n — 3, n — 2; i,j have the same parity}; 

2. Xi{n) and X 2 {n) have the same parity for every n > 0. 


In these terms we want to prove a lower bound for 

PhiMxiXin) = m*n). (A.15) 

Let Pi,qi,P 2 ,Q 2 denote the jump probabilities of X to the right, left, up and down, re¬ 
spectively. Notice that = pi — qi and m_ = P 2 — q 2 - 

Let H{n) denote the number of horizontal steps given by X in the first n steps. Then 
H (n) has a binomial distribution with success probability h := pi + qi. Let us assume that 
h < 1(2 for the remainder of the argument. A similar reasoning holds in the other case. 
Now given H{n) = k, we have that (A"i(n),X 2 (n)) = {Yi{k),Y 2 {n — k)), where Yi,Y 2 are 
independent simple random walks in Id with respective jump probabilities to the right 
= Pi/(Pi + Qi) and s = P 2 I{p 2 + 92 )- Notice that pi + qi > 0, P 2 + q 2 > 0. 
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We thus have that for any A: = 0,1,... , n 


^ PhiMA^A) = PA) = A 

= PhuhiA^^A) = m+APhiMA^An -k) = m-n)Ph^^h2APA) = A- (A.16) 

Let now Y((k) = (Yi(k) + k)/2. Then Y({k) is binomial with parameters k and r, and the 
hrst probability in (|A.16I) equals 


PhiMA^iA) = i'm+n + k)/2). (A.17) 

Similarly, the second probability in (IA.16I) equals 

PhiMA^iA -k) = {m-n + n- k)/2), (A.18) 

where Y^in — k) is binomial with parameters n — k and s. 

We will now choose k = kn either [hn\ or [hn\ +1 so that kn has the same parity as m+n. 
Notice that in this case kn > I and (m+n + kn)/‘2 is an integer. 

A straightforward recourse to Stirling shows that the last probability in ()A.16p is bounded 
from below by a constant times 1 / y/n. We next argue that the same holds for the proba¬ 
bilities in (IA.17P and (jA.lSh , and we will be done. We consider the first such probability; 
the second one can be similarly treated. 

We have that kn = {pi + qi)n + On-, where On G (—1,1], and (m+n -|- fc„)/2 = pin -|- 
Notice that Y((kn) has mean knV = pin -\- rOn- It again follows readily from Stirling that 
the probability in (IA.17P is bounded from below by a constant times 1 /An (notice that if 
Pi = 0, then On also vanishes). □ 


B Properties of the mean field free energy 


To prove Proposition [3] we study the free energy given by (I4.15P (two layers with a small 
vertical n.n. interaction). We shall exploit the smallness of e, observing that for e = 0 we 
have the well known explicit expression: 

/o(mi,m2) := - Ami) - ^ - ^(" 12 ), (B.l) 

where the entropy /(m) is given by 


Am) 

The function —\m? 
at 0 so that 


1 — m 1 — m l-|-m l-|-m 

■log—- —log——, mG[—1,1J. 


(B.2) 


2 ° 2 2 ° 2 
I (m) is a symmetric convex function of m with a quartic minimum 

/o(mi,m 2 ) > /o(0,0) -g c{m\ + m^). (B.3) 


Lemma 5. 

fAmi, m 2 ) - e < fAmi, m 2 ) < /o(m-i, m 2 ) + e. 


(B.4) 
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Proof. It follows at once from (lA.llI) and (jA.ll) since 7ro(-) — e < 7re(-) < 7ro(-) + e. 


□ 


Corollary 3. 


where 


There is c' > 0 so that 

inf A (mi, 777 . 2 ) 

mi ,7722 


Gc = I (mi, m2) e [- 1 ,1] 


= inf fe{mi,m 2 ) 

(B.5) 

(mi,m2)eC7^/ 

X [—1,1] : mj| < i = 1,2|. 

(B.6) 


Proof. Using (IB.4D and (|B.3p . we easily see that 

/^(mi,m 2 ) > A(0,0) +c(mf+ m|) - 2e, 

so that /e(mi,m 2 ) > /e(0,0) if (mi,m 2 ) ^ Gd with c' large enough, hence (IB.5p . □ 


We denote by /^(m), m = (m_|_,m_), the function /e(mi,m 2 ) when mi, m 2 are written in 
terms of m± as in (lA.Gp . In the sequel m = (m+,m_), h = (/7+,/7_) and 


{h, m) := /7+m+ + h-m^ = 2{mihi + m 2 h 2 ). (B.7) 

Given m and taking h as in (jA.Sp we then have 

fe{m) = -^{rn^ + rrif) + (/!>e(m) = -^(m, m) + ^(h,m) - Tre{h). (B.8) 


By ()B.5I) the inf of /e(m) is achieved in the set Gc for c large enough and the minimizers 
are critical points in such a set. Denoting by D the gradients, the critical points satisfy 

D/e(m) = —^m + D(/)e(m) = 0. (B.9) 

Then by Lemma [2] with 9 = m, 


m+ 


2e^ 


sinh(m+) 


cosh(m+) + e '^cosh(m_) 


m_ 


2 -e sinh(m_) 

cosh(m+) + e~^ cosh(m_) 


(B.IO) 


Of course m+ = m_ = 0 is a solution. In the next lemma we shall prove that any solution 
has m_ = 0. 


Lemma 6. For any x € M and any e > 0 the equation 


y 


2 e"" 


sinh(y) 

e^ cosh(x) + e~’^ cosh(?/) ’ 




has a unique solution: y = 0. 


(B.ll) 
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Proof. If y solves (IB.Ill) then so does —y. Therefore we only need to prove that there is 
no solution with ?/ > 0. Define 


U{y) := 


2 e- 


+ e ^ cosh(y) 


sinh(y). 


Since U{y) is not smaller than the right hand side of (jB.lip . the lemma will be proved 
once we show that U{y) < y for all y > 0. Suppose by contradiction that there is y > 0 
so that y <U{y). Then 

t . X _,sinh(y) 

e + e cosh(y) < 2 e -, 


which yields; 






2 n 


n>l 


n>l 


( 2 n + 1 )!' 


But this last inequality is not true, since + e ^ > 2e ^ and (2n + 1)! > 2(2n)! for 
n > 1 . □ 


We can thus put m_ = 0 in the first equation of (jB.lOjl which then becomes an equation 
for m_|_ alone. The proof of Proposition [3] is then a consequence of the following lemma: 


Lemma 7. There is 6 > 0 so that for all e > 0 small enough the equation 

sinh(x) 


X = 2e^ 


cosh(x) + e 


—, xG( 0 , 5 ) 


has a unique solution Xg and 


X, 


- V^\ < 


ce 


zji 


(B.12) 


(B.13) 


Proof. ()B.12p can be rewritten 


f , / N -F ,sinh(x) , 

e cosh(x) + e = 2e --—, xG(0, 5) 


X 


and therefore as F{x‘^) = 0 , where for z > 0 


and, 


vTl— 1 


n>2 


= 92{z) = 2e^Yl 


^n— 1 


n>2 


( 2 n + 1 )!' 


We have: 


(B. 14 ) 


F{z) = {e^ + e ^ + I + 2:51(2:)} - { 2 e^ + I e" + 2:52(2;)} (B. 15 ) 


F(12e) = {e^ + e-^ + 6ee^ + (12e)5i(12 e)}-{2 e^ + 4ee^ + (12e)52(12e)} 


= ce . 


(B.16) 
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On the other hand for any a < 1/6 there are (^o aiid eg positive so that for any 5 G (0,5o) 
and e G (0, eo) we have 

dF 

—— > a on (0, 6), 
dz 

Hence there is a unique z* G (0, (5) so that F{z*) = 0. Moreover 

ce2 

1^* _ i2e| < —, 
a 

proving the lemma. □ 


Theorem 6. For any e > 0 small enough the Hessian D‘^f^ of the free energy f^ is positive 
definite at the minimizers 


Proof. Since we shall be dealing with functions in the m and h domain, to prevent confu¬ 
sion we write and for the corresponding gradients (and for the Jacobian matrices). 
Analogously for the corresponding Hessian matrices. From ()B.9p and (lA.lip . by differen¬ 
tiating we have: D^f^ = —-I -|- and 

JO j. — ^ jH j. — 

Dm(Pe — 2 ’ ^raYe — ^ 

Since h{m) is the inverse of m[h): 


DhmDmh = I. 


On the other hand, by (lA.lOh . m = 2DhTT, so that 


D^m = 2G, G = D^tt 


and therefore 

2GDmh = I, Dmh = 

and, in conclusion, 

DlP, = ic-i 

Thus 

The elements of G are given in the equations which follow (IA.2I) and they must be com¬ 
puted at so that m_ = 0. Then G is diagonal and its diagonal elements, denoted by 
G++ and G _, are: 


2G++ 


2 gcosh(/i_|_)[e^ cosh(/i_|_) -|- e“^] — e*^ sinh^(/i+) 

cosh(/i+) -|- 

+ e~^ cosh{h+) ^ 2e^ e 

cosh(/i+)-|-e“^]2 “ cosh(h_|_)-|-e“^ “ 2' 


(B.17) 
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The last inequality holds for e small enough and it is proved as follows. We develop in 
Taylor series all terms up to first order in e, thus the equality below are meant modulo 
terms in e^. Recalling (jB.13p we also bound from below h\ = > 8e. The last fraction 

in (lB.17h is then bounded by 


2 + 2e ^2 + 2e_l + e^^ 

(1 + e)(l + hl/2} + l- e^2 + 4e“l + 2e- 

hence ()B.17D for e small enough. We have 


2G__ 


2e-" 

e*^ cosh(/i_|_) + e~'^ 



(B.18) 


for e small enough (and for any value of /i_|_). We have thus seen that I — 2G is diagonal 
and its diagonal elements are > | for e small. □ 


Corollary 4. For any e > 0 small enough there is c > so that for any C > 0 small 
enough: 


feim) - 


> cC^, for all m such that \m =F > (. 


(B.19) 


Proof. From what was already seen, the inequality holds if \m\ > with c large enough. 
The inhmum of fe{m) in \m =F rn^^'i\ > ( must then be achieved in the set 

{\m =F > C} Fl {\m\ < 


In such a set Dmfe / 0, thus the inhmum must be achieved at the boundaries, hence 

(1BT91) . 


□ 


C Multi-canonical constraints 

The setup is the following; I = Cq~ = [0,i-) nZ where, recalling (j2.2p . £_ = 
a > 0 and small. Let (mi(x),m 2 (x)) G [—1,1]^ for x G I, (mi(x),m 2 (x)) G [—1,1]^ for 
a; G Z \ /. Dropping the dependence on 7 and I, let 

T{m\m) = ^ 4 (mi(x),m 2 (x)) - ^ ^ J^{x,y)mi{x)mi{y) 

x£l 2 = 1,2 x^yGl 

+ X] J'rix,y)mi{x)mi{y)'^ (C.l) 

x^I 


where is the canonical free energy in Proposition [TJ 

Proposition [2] follows at once from the result below, which is the analogue for two layers 
of Theorem 6.4.1.1 of m. after the vertical interaction is added in. 
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Proposition 5. For all 7 small enough and all u = (ui, ^ 2 ), Ui G [—1,1], there is a unique 
rh such that rhi{x) = |/| Ui, z = 1 , 2 , and 

xGl 


F{m I m) > F{m \ m), for all m such that mi{x) = |/| ztj, i = 1, 2. (C.2) 

x£l 

Moreover the minimizer m is smooth in the sense that there is a constant c so that 

maxmax — uA < 07 “. 

i=l ,2 x&I 


Proof. The statement is trivially true when |ni| V |zz 2 | = 1, We therefore assume in the 
sequel that |zii| V |rt 2 | < 1. 

Remark. As in Section El for the free energy computation it is sometimes convenient to 
use the variables m±{x),x G I as in (IA. 6 P : we then write m(x) = (m+(x),m_(x)) and 
analogously m{x) = {m+{x),fh-{x)), and write 4>t.{mi{x), m 2 {x)) = (/>e(m(x)), with the 
function given by (IA. 8 I) . 

To find the minimizers under the above constraint, we introduce the Lagrange multipliers 
A = (A+, A_) and define 

Fx{m\m) = ^ (^(/>e(m(x)) - i(m(a:), A) - ^(K(x),m(x))^ 

X^I 

~ \ + (C.3) 

xj^yel 

with u± = U 2 ± ui, (a, b) = a+ 6 + + a-b- and 

i^±ix) :='^Jy{x,y)fh±{y), x£l. (C.4) 

y0 

Observe that F\{m\m) = F{m\rh) for all m under the constraint in ()C.2p . Let R = 
]Ti Ylxei We introduce an interpolating parameter s G [0,1] and define 

Fx,s{m\m) = ^{(/>,(m(x))-i(m(x),[A + ic])-s(^^{^([K(x)-K],m(x))} 

x^I xGl 

~ \ +^{\u)\i\ (C.5) 

xj^yel 


so that JA,i = d~\- To find the minimizer of | m) we need to find its critical points, 

namely the solutions of 


dFx,s{m\m) ^ ^ at \ o 

—BmPT “ =" 


where 


^{x) = + + k)J^{ x,y)m{y) 

yj^x,y&I 


(C. 6 ) 


(C.7) 
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By an explicit computation; 


am{x)chn[y) = ® 

which for 7 small enough is a positive symmetric operator, namely there is c > 0 so that 


sJ^(x, 2 /)(V'(x),V'(y)) > c^(V^(x),V^(x)) (C. 8 ) 

xdl xj^y£l xdl 


This shows that if there is a critical point of J-\^s it is unique and it minimizes By 

(|C. 6 p m is a critical point if D(j)t.{m{x)) = for all x. By LemmaO m{x) = 2D'Kf:{9{x)) 
which is (lA.lOp namely 


m+(x) 


2 e" 


sinh( 0 +(x)) 

cosh( 0 +(x)) + e~^ cosh( 0 _(x)) 


m_(x) 


2 e"" 


sinh( 0 _(x)) 

e‘^ cosh( 0 +(x)) + e~‘^ cosh( 0 _(x)) 


(C.9) 


Observing that the right hand side depends weakly on m as | Ylyei ^'ri^^y)^iy)\ — 
we then get, as in the proof of Theorem 6.4.1.1 in uni, that there is a unique solution mx^g 
of (jC.9jl which is the unique minimizer of J-a,s and whose fluctuations are of order 0 ( 7 “). 
To conclude the proof of the theorem it suffices to show that there exists a function A(s) 
such that 

= ^1^1’ for all s S [0,1]. (C.IO) 

xGl 

There is obviously a unique solution A(0) of (jC.lOl) when s = 0. We will find A(s) by 
solving the evolution equation 

= = 0> « e [0>1] (C-11) 

x£l x£l 


obtained by differentiating (jC.lOD and recalling that m = 2 Zl 7 rg( 0 ) =: 2'I'(0) where the 
explicit expression of 24' = 2(4'+(0),'I'_(0)) is given by the r.h.s. of (1C.91) . 

We will now prove that A(s) is differentiable and its derivative has order 7 “. We proceed 
by supposing that A is differentiable and get a formula for its derivative. We will then 
check that the primitive for such an expression is indeed A(s). The formula will also show 
that A'(s) has order 7 “. 

By (|C.lip we first need to show that 0x[s),s is differentiable in s. But 9x(s),sA) defined 
through (IC.7P can be expressed as 

= + (C.12) 


where 

A(s) 


We have: 


i(A + «;) and 


= i (k(x) -k)+ 

|0i(l)(x)| < C7“. 


x,y&I 


Jy{x,y)mx{s){y) 


so that 
(C.13) 
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with: 


(C.14) 


d(^\(s),s 

ds 


ds 

dxOx{s),s{^) 


(x) = dxOx{s),s{^)j^ + 

(x) 


dX dOxis),. 


ds ds 


■(x) 




dxO^'^'> + s dxO^^'^ = + s ^ J^(x, y)Dxmx, 

yGl 


(C.15) 

(C.16) 


where dx ■= 


'a := ® identity matrix, and mx is defined by equation: 

m{x) = ^{ex,m{x)) (C.17) 

where 9x,m as the same expression of 9 (jC.7p but as a function of A, m, namely: 


~ 1 s 

Ox,m{x) = -{X + R)-- 


{k{x)-k)+ ^ Jy{x,y)m{y) 

yjtx,yel 


(C.18) 


We postpone the proof that mx is differentiable in A and there is a constant c so that: 

|iAAm| < cy" (C.19) 

which, recalling (|C.16p . implies that there is a constant cj so that: 

l-cjr<\\dx9\\<l + cjr- (C. 20 ) 


Going back to (IC.lip and denoting by the operator that acts on a vector v{x) as 
Tgv{x) := Ylxei ^9'^(dx{s),s{x))v{x), we write in a compact form: 


,dX 


There are positive constants cf (see ()C.9D 1 so that: 

\I\c- < ||T,|| < |/|c^ 


(C.21) 


(C.22) 


Then, by ()C.20p . [T6I(9 a6*a(s),s] is invertible if 7 is small enough. Hence ^ is well defined 
and 


dX 


ds 


Last inequality follows by (1C.151) . 


'^edx9x{s). 


1 -1 


Te9 


( 1 ) 

X{s),s 


< C-f 


(C.23) 


Let us now prove (IC.19p . Differentiating ()C.17I) : 


Dxm = DgT 


dx9 + dyn9 ■ Dxm 


(C.24) 
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then 


1 - • ( 9^6 


D\m = Dg^ ■ d\6 


(C.25) 


By (IC.ISP we see that d\6 = 1/21, and D„iO is a smooth function order 0 ( 57 "). Then for 
7 small enough WDg'^DmOW < 1 and 1 — Dg^DmO 


is invertible and 


D\m = 


1 - Dg^DmO 


-1 


Dg^ • d^e 


is well defined and there is a constant c so that: 

\D\m\ < C 7 ° 


(C.26) 


(C.27) 


D A contraction property of the mean field free energy 


In this section we shall prove Proposition 01 The basic bound comes from the analysis of 
the previous section. As we use here the variables (mi,m 2 ) instead of (m+,m_) we need 
to translate the results into the new variables. The minimizer = (m^^, 0) becomes 
(me,me), rrie = j^. We have 


Lemma 8 . Let G++ and G _ be given by (IB.lTp - dB.lSp . then 


52 




(me,md = G++ + G —, i = 1,2, 






-(m£,md = G++ - G—. 


Both G++ + G _ and G++ — G _ are non negative, and for i = 1,2 

M I , M ^ f 

+ 1^ - 2- 


(D.l) 


(D.2) 


Proof. 


dir^ StTc StTc drr^ ^ 


dhi dh^ dh- dh2 (9/i_ 

|2.^ 02.^ a2.^ q2 


+ -mf- - 2 




92 vr,; 9271^ 


9/i2 9/i^ 9/i?_ dh+dh- ’ dh2dhi 9h+ 9/i?_ 

hence ()D.ip because (2mg,0) = 0. 


□ 


It follows by continuity that: 
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Corollary 5. There are cq > 0 and r < I so that the following holds. Call 


Ri,j = sup 

h=(hi me |<co,^=l,2 


dhidhj 


{h)\. 


Then 

Ri,j <r, i = l,2 

j=l,2 

and the matrix {1 — R) is invertible. 


(D.3) 


(D.4) 


We are now ready for the proof of Proposition [H We fix throughout the sequel a pair 
(x, i) and (x, i') of vertically interacting sites in Ain, thus (x, i') = Vx^i and use the above 
properties to study the function g^{m) introduced in (|5.4I) . We write m = {mi,mii) and 
write a, = ax,i, aii = Oxy, the latter defined in (14.Sp . We also write A* and Aj/ dropping 
the superscript u on which they depend via ()5.5p . We finally shorthand aj for Oxj. 

Lemma 9. There is a constant c so that for any u G Rfx,i,i', see (I5.3p . 

^-^<C + c7“, j = (D.5) 

l-Oj 

where A®'^ is defined in Proposition 
Proof. By (I5.5h 


^_ j 

y¥^^-{yd)&C^ 

+ X] Ji{x,y)iu{y,j) -mf). 

y-[yd)t{ci~ ’■’uAo} 

We add and subtract J^{x,y) to J-y(x,y) where J^{x,y) is obtained by averaging J.y{x,y') 
over Cy~’R In the term with J.y{x,y) we can replace u{y,j) by its average and use (|5.3I) . 
(j4.8p and that is a probability kernel to get the bound (1 — aj)Q. The sum over the 
terms with J.y{x,y) — Jj{x,y) is bounded by c' 7 ", by the smoothness of Jy. 

□ 


In the sequel we shall only use the bound (jP.Sp and not the specific form of the Xj. 
By differentiating (j5.4jl we get 


d'^ge 

dmjdmji 


dmjdniji ’ 


j,j' G {i,i'}. 


^ ~ dm dm., diagonal with entries — aj, j = i,i'. The minimum of —I" is 

at 0 and —I"{R) = 1 . Since Oj <1/2 (this follows from the choice of x and the symmetry 
of J; see (14.8p ). we then conclude that: 


Lemma 10. There is ci > 0 so that g^ is strictly convex for e < ci and for any such e it 
has a unique minimizer m (called in Proposition^ . 
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In the sequel we tacitly suppose e < ci. The critical point m of satisfies 


D^^{m) = 6 


(jiiUii T T Aj/) 


(D. 6 ) 


Then, by Lemma O 


m = DhTTeiO), 


(D.7) 


where Dh^e{(^) is the gradient of 'ne{hi-, / 12 ) computed at (hi, /i 2 ) = 0. We shall study (ID.7P 
distinguishing among the possible values of a, and Oj/ which depend on the horizontal 
distance of {x,i) and respectively (x,i') from Aq. Because of the geometric properties of 
Ain, only three cases can occur: (i) Oj = Oj' = 0; (ii) Oj = Oj' G (0, |]; (iii) a* G (0, 
and di' = 0 or viceversa. Case (i) occurs when the horizontal distances of (x, i) and (x, i') 
from Aq are both > 7 “^. Case (ii) is when the distances of (x,i) and (x,i') from Aq are 
both < 7 “^ and case(iii) is when one is < 7 “^ and the other > 7 “^. We start from case 
(i) which is the easiest. 


Case (i). By (ID.SP for 7 small enough \0j — nie 


< Co, Co as in Corollary [5l hence by (ID.SP 


rui — nif 




jj — 


< r < 1 


in agreement with in Proposition [4] after setting j') = Rjj' and recalling that in 

case (i) dj = 0, j = i,i'. 

Case (ii). 9 in (jD.7p is now (after adding and subtracting m^) 

9j = rUe + d{mj — m^) + Xj — (1 — d)m^, j = i, i' . (D.8) 

Since 6 depends on m, (ID.7P is an equation in m and not a formula for m as in case (i). 
We introduce an interpolating parameter t G [0,1] and define 


Ojit) 


+ d{mj 


me) +1 




j = l,2 


(D.9) 


calling m{t) the solution of (jD.7jl with 6 replaced by 9{t). Observe that m(0) = me is the 
solution at t = 0 while the solution at t = 1 is what we want to find because 6{1) = 9. 

Supposing that m{t) and its derivative m{t) exist we can then differentiate (ID.7P to get 


mi = 


E 

p=i,i' 


Kjpidrhp + (Ap - (1 - d)me)}, 


K 

dhidh. 


■m) 


(D.IO) 


where 9{t) is computed at m = m{t). If moreover \m{t) — me\ < 2(^ then \9j{t) — me\ < 2<^ 
and I — dK is invertible and we have 


m = V{m,t) := {1 — dK) ^K{X — {1 — d)u^'^), = (m^, m^), A = (Aj, A*/). (D.ll) 

The evolution equation (ID.lip starting from m(0) = me has a unique solution till the 
hrst time T when |mj(r) — me\ = cq, because by Corollary[5l (1 — K{t)) is invertible and 
smooth for t <T and we have 

00 

\mj{t) — me\ < t ^ a" ^ {R"'~^^)jp\Xp — (1 — d)me\. (D-12) 

n=0 p=i,i' 
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Set 


OO 




3,3 ■ 


n=0 


Then, by ()D.4jl . we get 


Cx,i,i'{j,j') < (1 - a)- 


— ar 


<r, j = i,i 


Thus \mj{t) — rriel — for t < min{T, 1}, hence the above holds till t = 1 and Proposition 
mis proved in case (11). 

Case (111) with a* = a > 0 and a*/ = 0 (same proof applies when a* = 0 and Oj' > 0). Here 
9 = + a(mj — m^) + A, — (1 — a)m^, rrif: + (A*' — 

and proceeding as In case (11) we set 

9{t) = (jrif: + a{mi — me) + t[Xi — (1 — a)me],me + t(Aj/ — rrie)^ . (D.13) 

Analogously to (|D.10h . 

rhi = Ki^iarhi + {Ki^i{Xi - (1 - a)me) + Ki^ii{Xii - me)}. 

Since Ki^iU < 1 till when |mi(t) — me\ < 2( proceeding as in case (11) we get that the 
evolution equation has solution till time t = 1 and 

OO 

|m,:(l) - me\ < *)""{/?*,i|Ai - (1 - a)me\ + Ri,i'\Xi> - me\}. (D.14) 


n=0 


We then set: 


Cx,i,Ahi) = Cx,^,Ai,i') = 


Ri 


1 — aR 


' 2,2 


1 — aRj 


which verifies the condition in Proposition [4] because 

E , , R^ j (1 — Qj) Ri V T — ClRi i 

Cx i i) = - - H- - — < - — < r. 

l-aRi^i 1 - aRi^i - 1 - aRi^i 


(D.15) 


Since rui/ = ^^{0), 

Olljf 


|mj/ — me\ < RiiAA' — 'nie\ + Ri',i{a\mi — me\ + |Ai — (1 — a)me'^ ■ 


We then set 


Ri 


Cxjiji'ii ) O Ri',i' T 0,Ri',i „ 

i 

n i ■! — T} ( ~ A I 1 

Cx i i' (^ ) *) Ri' ,i Tj Tl 0, 

V i — aKi i 


(D.16) 
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and 


. -I- 


- a) 


+ (1 — a 


1 (xRi^i 

Ri , i ' -^2,2(1 


1 — 1 — 


< r — aRii^ii 1 — 


Ri , i ' “ 1 “ -^2,2(1 


1 (iRi^i 

having bounded in the last bracket Ri^i/ + Ri^i < 1. 


< r 
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